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INTRODUCTION 

The thermodynamic calculations of the decarburiza¬ 
tion of liquid metal, in particular, the estimation of oxy¬ 
gen distribution for the oxidation of carbon, iron, and 
other melt components, assume the determination of 
the carbon and oxygen activities in a metallic solution 
at given composition and temperature. 

For dilute binary solutions of carbon in liquid iron, 
zero-order (lny c ) and first-order (£ c ) molar self-inter- 
action parameters are used. In the case of a relatively 
high carbon concentration, a second-order (p c ) inter¬ 
action parameter should also be taken into account. The 
values of the parameters available in reference books 
[1-3] (published in different years) differ substantially. 
This makes the use of the data for the calculations dif¬ 
ficult. Thus, to increase the accuracy of the calcula¬ 
tions, the correct selection of these magnitudes and the 
determination of their temperature dependences are 
necessary. 

When performing thermodynamic calculations of 
the carbon activity in alloyed melts, equilibrium oxy¬ 
gen contents should be determined and taken into 
account along with the concentrations of other compo¬ 
nents. 


where A Fe and A c are the atomic weights of solvent 
(iron) and solute (carbon), respectively. 

According to [1], the Gibbs energy of reaction (1) is 

AGc = 22600 -43.27\ J/mol. (3) 

Equations (2) and (3) can be used to determine the 
temperature dependence and activity coefficient of car¬ 
bon in its infinitely dilute solution in liquid iron: 

InYc = 2.013, Yc , 1873 k = 0.57. (4) 

In later work [2], the temperature dependence of the 
Gibbs energy of reaction (1) was refined: 

AGc = 17230 - 39.87 T, J/mol. 

This allows us to refine the temperature dependence 
of Yc : 

oo 2072 

l n Yc = —jr- ~ 1-726, Yc.isvbk = 0.54. (5) 

The values of y F determined by Eqs. (4) and (5) dif¬ 
fer slightly; nevertheless, Eq. (5) is preferable. 


ACTIVITY COEFFICIENT OF CARBON 
IN AN INFINITELY DILUTE SOLUTION 
IN LIQUID IRON 

The Gibbs energy AG( of the reaction of transition 
of carbon from the standard state “pure substance” into 
the standard state “hypothetical 1% solution” 

C(graphite) ~ [C]i% (1) 

allows the activity coefficient of carbon in an infinitely 
dilute solution y( : to be determined by the equation [4] 


FIRST-ORDER SELF-INTERACTION 
PARAMETERS OF CARBON IN LIQUID IRON 

The temperature dependence of the first-order self- 

c 

interaction mass parameter of carbon in liquid iron e c 
and its value at 1873 K are given in [1]: 

ec = ^ + 0.0581, 4 m3 k = 0.143. 

The conversion of the mass parameter e c into the 
molar parameter £ c by the formula 


RT lnyc = AG^-RT 


^Fe 

100A C ’ 


( 2 ) 



c 
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+ 



( 6 ) 
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Fig. 1 . Fe-C phase diagram [5] and the temperature dependence of the parameter p c calculated from these two versions of the 
phase diagram. 


gives e c = 7.8. However, the authors of [1] in [3] rec¬ 
ommend using the highest value of this parameter, 

namely, £c,i 873 K = n - 

The conversion into the mass parameter by the for¬ 
mula 

C _ 1 A F c A c - A Fe 

c 230 A c c A c 

gives ^c, 1873 k - 0.21. 

The refined value of the e c parameter at 1873 K, 

e C, 1873 k = 0.243, 

is given in [2]; however, its temperature dependence is 
not reported. 

The temperature dependence of the e c parameter 
given in [1] can be used; however, its value at 1873 K 
should be increased from 0.143 to 0.243. In this case, 


the temperature dependence of the e c parameter takes 
the form 


Using this dependence in Eq. (6), which is applied 
to calculate the molar parameter from the mass param¬ 
eter, we obtain the temperature dependence of the 
molar parameter of carbon self-interaction: 

c 7787 c 

e c = ~y~ + 8.622, e£ I873 K = 12.8. (7) 

ESTIMATION OF SECOND-ORDER 

SELF-INTERACTION PARAMETERS 
OF CARBON IN LIQUID IRON 

Let us estimate the second-order interaction param- 
eters of carbon p c using two versions of the Fe-C 
phase diagram (Figs, la, lb) [5]. We assume that, for 
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c 

Table 1. Results of the calculation of the parameter p c from the Fe-C phase diagram 


T, K 

From the diagram in Fig. la 

From the diagram in Fig. lb 

c -19733 . 

p c - T +1.33 

[C] sat . % 

, Scit 

c 

Pc 

x c , sat 

c 

Pc 

1426 

4.25 

0.171 

-12.62 

0.17 

-12.08 

-12.51 

1473 

4.35 

0.175 

-12.14 

0.175 

-12.14 

-12.07 

1573 

4.57 

0.182 

-10.8 

0.184 

-11.69 

-11.21 

1673 

4.79 

0.190 

-10.37 

0.191 

-10.77 

-10.46 

1773 

5.00 

0.197 

-9.84 

0.197 

-9.84 

-9.8 

1873 

5.20 

0.203 

-9.32 

0.203 

-9.32 

-9.21 

1973 

5.37 

0.209 

-8.83 

0.207 

-8.19 

-8.67 


the saturated solution of carbon in liquid iron, the car¬ 
bon activity at the liquidus is a c sat = 1. In this case, 

Inflc.sat = ln* C ,sat+ lny c ,s at = 

where 

lny C ,sat = l n Yc + £c*C, S at+ Pc(*C, S at) 2 - 

From these equations, we can derive an equation for 

Q 

the calculation of the p c parameter: 

c -(In Yc + £c*c, sat + ln*c, sat ) 

Pc = --■ (8) 

V-^C, sat/ 

The parameters calculated for a temperature of 
1873 K by Eqs. (5) and (7) are lny^ = -0.620 and = 
12.8, respectively. According to the phase diagram in 
Fig. la, the point at T- 1873 K corresponds to the com¬ 
position [C] = 5.20% or x c = 0.203. In this case, the 

parameter calculated by Eq. (8) is p c = -9.32. 

The analogous calculations of the p c parameter 
were performed at other temperatures of the melt, i.e., 
at 1426-1973 K, using points at the liquidus lines in 
both versions of the Fe-C phase diagram (Figs, la and 
lb). Table 1 gives the carbon concentrations corre¬ 
sponding to these points and the calculated values of 

the p c parameter. 

According to the dilute solution theory in [6], the 
temperature dependence of the interaction parameter 

should take the form p- = ^ + b. The p£ parameters 

given in Table 1 adequately fall in a segment of the 
straight line plotted on the p c - 1 IT coordinates 
(Fig. lc). The least squares processing of the data in 
these coordinates leads to the regression equation 

c -19733 

Pc(T) = -i^ + 1.33 (9) 


with a high correlation coefficient (r - -0.981). 

Q 

The values of the p c parameter calculated by 
Eq. (9) (Table 1, the right column) are comparable with 
those calculated from both versions of the Fe-C phase 
diagram. We suggest using Eq. (9) for the estimation of 

the p c parameter in a temperature range of 1426-1973 K. 

Table 2 lists the self-interaction parameters of car¬ 
bon in liquid iron and their temperature dependences 
calculated using reported data and the Fe-C phase dia¬ 
gram in [5]. 

With the parameters given in Table 2, we can calcu¬ 
late the activity coefficient and activity of carbon in liq¬ 
uid iron at a given concentration: 

lny c — lny c + £c x c + Pc( x c) > a c ~ Yc x c- 

At a low carbon concentration, the second-order 
parameter p c may be neglected, and the activity coef¬ 
ficient and activity of carbon can be calculated by the 
equations 

lnyc = lnyc + £c*c; a' c = i c x c . 


Table 2. Self-interaction parameters of carbon in liquid iron 


Value at 1873 K 

Temperature dependence 
of the parameter 


, <*> 2072 , 

OO 

cn 

in 

o 

II 

1 u 

lny c - -y--1.726 

c 


8c = 12.8 

—-h 8.622 

c T 

c 

c -19733 , 

Pc =-9.21 

Pc- T +1-33 
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Fig. 2. Relative error of calculating the carbon activity in 
liquid iron at T = 1873 K that is related to the neglect of the 
C 

p c parameter. 


five boiling depth h ~ 15 cm. When calculating, we 
took into account the effect of alloying components 
(chromium, nickel), the residual silicon and manganese 
contents, and the equilibrium oxygen content. 

We used the reaction of carbon oxidation at 1500- 
2000 K, 

^-■(graphite) ^ ^“^2(g) — CO(g), 

AG° = -114 515-86.0157; 

mol 

whose equilibrium constant is 

AG° 

„ P CO 0.1 RT 

K - - - - - e 

a c a o x c x oYcYo 

At given temperature and carbon concentration x c , 
the carbon activity and equilibrium oxygen concentra¬ 
tion are 


a c ~ Yc x c> x o 


0.1 

^" x oYcYo 


( 10 ) 


The related relative error of the activity calcula¬ 
tion is 


The activity coefficients of carbon y c and oxygen y 0 
are expressed in terms of interaction parameters £■: 



it increases with the carbon concentration (Fig. 2). 

At a carbon concentration [C] = 0.1% in the melt at 
1873 K, the relative calculation error is negligible (8 C = 
0.02%); at [C] = 0.1%, the error is 5 C = 1.88%. In the 
saturated solution ([C] = 5.2%), the error reaches 45%. 

When calculating the carbon activity in liquid iron 
with a carbon content higher then 1%, it is reasonable 

to use the parameter pc . 


CARBON AND OXYGEN ACTIVITIES 
IN AN ALLOYED STEEL 

At present, an innovative (vacuum) technology for 
deep decarburization of alloyed melts in complex out- 
of-furnace treatment plants is being implemented at 
high-quality metallurgical works [7]. In connection 
with the use of such innovative technologies, it is a 
challenge to analyze changes in the carbon and oxygen 
activities, from the onset of the oxidation process at a 
sufficiently high carbon concentration (0.40-0.70%), to 
the end of the process, when the weight fraction of car¬ 
bon decreases to several hundredths of a percent. We 
calculated the activity coefficients and activities of car¬ 
bon and oxygen for five pilot heats of a corrosion-resis¬ 
tant steel, which differed in metal composition and the 
temperature at the beginning and end of an oxidation 
process in a vacuum (Table 3). The partial pressure of 
CO in gas bubbles during metal boiling was taken to be 
Pco = 0.1. The principal contribution to this value was 
made by the ferrostatic pressure of the metal at an effec- 


lny c = Inyc + e c x c + e c x o + E 


Cr 

C X Cr 


Ni 

+ £r X\ 


Mn 


lny 0 = Inyo + £o x o + e o x c + £o x cr 

Ni Mn Si 

+ £ 0 x Ni + £ 0 x Mn + £ 0 x si, 


( 11 ) 

( 12 ) 


where 


Inyc 


Inyo = 


2072 

T 

14424 

T 


1.726, 

+ 3.275, 


Yc, 1873 K _ 0.538, 
Yo, 1873 K = 0.012, 




7787 

T 


+ 8.622, 


o 

£o 


115000 

T 


+ 50.66. 


The interaction parameters at T = 1873 K [2] are as 


follows: 

C 

O 

Cr 

Ni 

Mn 

Si 

i 

e c 

12.8 

-20 

-2.0 

2.37 

-1.88 

9.75 

£o 

-20 

-10.7 

-4.84 

2.40 

-4.7 

-7.10 


At other temperatures, the parameters of interaction 
of carbon and oxygen with alloying components were 
corrected by the equation 

i _ j 1873 

£ l, T ~ £ i, 1873 K j ■ 


The calculation, according to the system of 
Eqs. (10)—(12), is difficult owing to the fact that, to cal- 
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Table 3. Compositions of alloyed melts upon vacuum decarburization in a 10-t unit and the activity coefficients of carbon and oxygen calculated for these compositions 


INTERACTION PARAMETERS OF CARBON IN LIQUID IRON 


549 



RUSSIAN METALLURGY (METALLY) Vol. 2007 No. 7 


* UM means unmeasured. 
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Fig. 3. (I) Dependence of the activity coefficient of carbon 
y c on its molar fraction Xq for the averaged composition of 
an alloyed melt and (2) the equilibrium curve for an Fe-O- 
C melt. Points in the curve correspond to the activity coef¬ 
ficients of carbon y c in the alloyed melts of each heat calcu¬ 
lated from the metal composition and temperature. 


«C 



Fig. 4. Dependence of the activity of carbon on its molar 
fraction for the averaged compositions of (7) alloyed and 
(2) Fe-O-C melts. Points in the curve correspond to the 
activity of carbon in the alloyed melts of each heat calcu¬ 
lated from the metal composition and temperature. 


To 



Fig. 5. Dependence of the activity coefficient of oxygen Yo 
on the molar fraction for carbon for the averaged composi¬ 
tions of ( 1 ) alloyed and (2) Fe-O-C melts. Points in the 
curve correspond to the activity coefficient of oxygen in the 
alloyed melts of each heat calculated from the metal com¬ 
position and temperature. 


x c x 0 , 10 



Fig. 6. ( 1 ) Dependence of the product of the molar fractions 
of oxygen and carbon on the carbon molar fraction for the 
averaged composition of an alloyed melt and (2) the equi¬ 
librium curve for an Fe-O-C melt. Points in the curve cor¬ 
respond to the product of the molar fractions of oxygen and 
carbon in the alloyed melts of each heat calculated from the 
metal composition and temperature. 


culate a c and x Q by Eq. (10), x Q in Eqs. (11) and (12) 
must be available. That is why we used the method of 
successive approximations. 

The activity coefficients and activities of carbon and 
oxygen were calculated using the alloying component 
contents in the melt averaged over five heats (17.41% 
Cr, 6.355% Ni, and 0.89% Mn) and the alloying com¬ 
ponent contents in the Fe-C-0 melt at an average tem¬ 
perature of 1945 K. The calculations were performed 
using the Microsoft Excel 2002 software package. 


The initial data, the calculated activity coefficients 
and activities of carbon, and the equilibrium values of 
y 0 , x q , and a Q for oxygen are given in Table 3. 
Figures 3, 4, 5, and 6 illustrate the calculation results. 
The calculation of the equilibrium curves of the quanti¬ 
ties for the average composition was performed to 
achieve the minimum molar fraction of carbon, which 
corresponds to the maximum equilibrium solubility of 
oxygen. For the alloyed melt, x c min = 0.00093; for the 
Fe-C-0 solution, x c min = 0.000049. 
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As is seen from Fig. 6, the product of the molar frac¬ 
tions of carbon and oxygen x c ■ x Q in the alloyed melt is 
located substantially higher than that for the iron-car¬ 
bon melt. The existence of a minimum in the product of 
the molar fractions of carbon and oxygen, which was 
found by Filippov on experimentally studied iron-car¬ 
bon melts [8], is confirmed by the thermodynamic cal¬ 
culations. In the presence of alloying elements, the 
value of the minimum increases and it shifts toward 
high carbon concentrations. 

The low activity coefficients y c and y 0 i n alloyed 
melts as compared to those in iron-carbon melts lead to 
a decrease in the activities of carbon and oxygen in the 
alloyed melts. 

Alloying components such as chromium and man¬ 
ganese are shown to hinder deep decarburization of 
alloyed melts. However, the use of vacuum during 
decarburization allows a carbon concentration of 0.02- 
0.04% to be reached in melts containing 17-18% Cr at 
the end of the oxidation stage. 
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